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LUNAR MISSION PERFORMANCE EVALUATION PROCEDURE FOR 
ORBIT -LAUNCHED NUCLEAR VEHICLES 


By Robert E. Austin and George B. Kearns 


SUMMARY 


An estimation procedure, based primarily on the analysis of 
Reference 1, is described and reference data are presented by which the 
influence of certain design parameters and lunar mission velocity 
requirements can be determined for orbit -launched nuclear vehicles. 

The use of this procedure allows a rapid, unsophisticated method of 
calculating the performance of nuclear propulsion systems for specific luna 
mission profiles. The performance parameters considered are as 
follow: payload ratio, characteristic velocity requirement, propellant 

ratio, flight path angle, central angle, burning time, and injection 
altitude. The effects of initial thrust -to -we ight ratio, specific impulse, 
lunar transfer time, and initial orbital altitude are considered. The 
range of specific impulse assumed is restricted to that considered 
practical for nuclear heat exchanger propulsion systems. 

Application of this estimation procedure is demonstrated by 
numerical examples included in Appendix A. The results of these 
and other examples are compared with exact computer solutions and 
are shown to be well within the accuracy usually required in preliminary 
design studies. 



SECTION I. INTRODUCTION 


With increasing interest in the application of nuclear rocket 
space systems to lunar transportation, it becomes important to have 
available a rapid, practical method for relating nuclear vehicle design 
and performance parameters to trajectory characteristics, once 
preliminary mission profiles and conceptual vehicle designs are 
established. A comprehensive procedure of this nature can provide a 
useful tool for conducting feasibility studies of contemplated missions. 

The estimation procedure used must have sufficient accuracy for 
preliminary design studies. 

The purpose of this report is to present reference data and a 
procedure by which the effect of certain design parameters, upon the 
basic performance of orbit-launched nuclear vehicles, can be 
determined. Lunar transportation systems will require large initial 
gross weights in earth orbit. Conventional chemical vehicles are 
comparatively inefficient for the orbit-launch boost phase of such missions. 
Hence, the range of specific impulse assumed in this report is restricted 
to that considered practical for nuclear heat-exchanger propulsion systems, 
namely 700 to 1000 seconds. While the reference data presented are 
restricted to nuclear rocket systems, the general method of performance 
estimation can be applied with equivalent success to any orbit-launched 
system. 


The method of calculation of the parametric data presented in this 
report is primarily based on the analysis of Reference 2; and, the performance 
evaluation procedure is based on that of Reference 1. Since the scope of 
Reference 1 is limited to escape and planetary missions, this report extends 
that work to include a range of lunar missions with transfer times of 
49- 75 hours (earth escape) through 120 hours. Although the work done in 
Reference 2 covers lunar missions, it is concerned primarily with the 
influence of structural weight parameters on the optimum thrust -to-weight 
ratio of orbit-launched vehicles. Hence, the method of calculation is 
briefly presented in this report for ease of reference. 

The reference transfer time assumed in this investigation is 
70 hours, because it is representative of current lunar mission studies. 

Only transfer times between 49. 75 hours (earth escape) and 120 hours 
have been considered. Performance evaluation procedures for transfer 
times less than 49- 75 hours can be found in Reference 1, where the 
missions correspond to earth escape and hyperbolic excess speeds. 
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It is pointed out in Reference 2 that tangential thrust vector 
orientation is near optimum for near escape, escape, and particularly 
for hyperbolic excess speeds; hence, tangential thrust is the only mode 
considered in this report. A single stage with constant thrust is assumed 
from earth orbit to final injection. Initial thrust -to -weight ratio, F/W 0 , 
is varied from 0. 09 to 1. 0, with both limits being arbitrarily selected. 
Initial orbit altitudes are assumed to lie between 185. 2 and 740. 8 km 
(100 and 400 N. M. ). 

Reference curves for each of the parameters are presented for 
the specific impulse range considered, along with correction curves for 
variations in transfer time and initial orbit altitude. Variations in vehicle 
structural and engine weights can be accounted for in one algebraic 
equation. 

Calculations were made by numerical integration on the IBM 7094 
digital computer. The equations of motion are related to a spherical 
earth. The analysis is limited to planar trajectories with all aerodynamic 
and perturbative forces being neglected. 


SECTION II. ASSUMPTIONS 

The basic assumptions made in the analysis are summarized 
as follows: 

1. Acceleration of a single stage out of a circular orbit, with 
constant tangential thrust. 

2. Range' of independent performance parameters: 

I = 700 - 1000 sec 
sp 

T = 49. 75 - 120 hr 
tr 

h = 185.2 - 740. 8 km (100 - 400 N. M. ) 
o 

F/W q = 0. 09 - 1. 00 
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3. Reference conditions: 

T = 70 hr 
tr 

h = 555. 6 km (300 N. M. ) 
o 

A = B = 0. 1 
C = 0 

4. Mean spherical earth model with: 

km 3 

^= 398,613 ^ 2 . 

r = 6371. 2 km (3440. 2 N. M. ) 

© 


SECTION in. ANALYSIS 


The impulsive velocity requirements, for injection into a 
transfer orbit having 80 -km periselenum, were computed using 
Figure 1. The velocities in Figure 1 were calculated assuming 
horizontal injection (6 = 90°), whereas the flight path angle at injection 
in this analysis was not restricted to this condition. It is pointed out 
in Reference 2 that the errors introduced by this approach are negligible. 
The method of calculation is based on the sphere -of- influence concept, 
reducing the n-body problem to two 2 -body problems. 

Since the thrust is assumed to be tangential (j3 =' 0°), no attitude 
restriction can be placed on the resulting trajectory. Hence, the 
attainment of the desired transfer time injection velocity is assumed to 
be a sufficient end condition. 


From Reference 1, the general equations of motion for a vehicle 
leaving circular orbit are: 


_ F cos /3 
m 

A _ F sin g 
mV 


- g cos 9 
+ (y - - T ) Sin 6 


( 1 ) 

(2) 


4 



Numerical integration of the equations of motion determine the 
velocity and flight path angle. The range and altitude are then calculated 
by: 

r 

x = f — V sin 6 dt 
J r 


h = h + f V cos 0 dt 
o J 

The central angle is determined by: 


+ =/ 


dt 


© 


It should be noted that the initial thrust -to -we ight ratio is hidden in 
Equation 1, since 


F 

m 


= ( 


W 


) g 


n 


In Reference 1 it is pointed out that the weight characteristics of 
a vehicle must be defined before an evaluation of payload performance 
can be made. Defining the weight characteristics, as in Reference 1, 
the first group, A, is composed of the engine, propulsion system 
hardware, and structural members which may be assumed proportional 
to thrust. The second group, B, is composed of propellant tankage and 
any propellant residuals which may be assumed proportional to 
propellant loading; and the third group, C, contains interstage structure, 
astrionic gear, and various miscellaneous equipment which may be 
assumed proportional to the initial weight of the vehicle. The net 
effective structural weight is expressed as 

T = AF + BW + CW 
N p o 

Dividing the above equation by W Q to nondimensionalize, gives 
e N = A (F/W o ) + B& + C (3) 

The relation between stage payload, propellant loading and net 
structural weight is 

W T = W - W - T 
L o p N 
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Nondime ns ionalizing, 


X = 


- C - e. 


N 


( 4 ) 


Rearranging Equations 3 and 4, 

X = 1 - A (F/W q ) - (B + 1) C - C (5) 

where: X - W fW 
L o 


i = W /W 
P o 


The propellant ratio is defined for constant thrust and specific 
impulse by 


c - 


_F 

W 



( 6 ) 


All reference curves shown for payload are for the specific 
weight parameter C equal to zero. The effects of non-zero values of C 
are shown to be accounted for by a single algebraic equation later in 
this report. 

Reference values for initial orbit altitude, transfer time, and 
specific weight parameters were assumed arbitrarily, but should be 
reasonable since they only provide a basis for the estimation procedure. 
The reference values chosen for this report are: 

1. T = 70 hr 
tr 

Z. h = 555.6 km (300 N. M. ) 
o 

3. A = B = 0. 1 

4. C = 0 


Reference curves are shown in Figures Z through 8 for the various 
parameters. For preliminary design purposes, payload and characteristic 
velocity are perhaps the most important of all the performance parameters. 
Hence, the curves of payload ratio and characteristic velocity shown in 
Figure Z and 3, respectively, are plotted on a highly expanded scale. 
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Once reference data are delineated for the individual parameters, 
consideration must be given to variations in vehicle specific weights, 
transfer time, and initial orbit altitude. The effects of such variations 
and the means of estimating each of the parameters will now be discussed. 

A. PAYLOAD RATIO 

Of all the parameters considered, only X is affected by 
changes i.n A and B. For reference conditions, Equation 5 can be 
written 


X 

re 


f 


- A (F/W ) - (B + 
ref ° ref 


1) c 


ref 


(7) 


where both V , and L are functions of F/W„ and I , as illustrated 
ref ref u sp 

in Figures Z and 4, respectively. Arbitrary values of A and B can be 

expressed in the form, 

A = A + AA 
ref 

= 0.14- A A (8) 


B = B _ + AB 
ref 

= 0. 1 + AB 


(9) 


Hence, for arbitrary values of A and B, Equation 7 becomes 




ref^A, B 


1 - A (F/W q ) - (B + 1) ^ ef 


where the terminology is adopted that subscripts outside parentheses, 

or brackets, enclosing one of the parameters indicate non-reference conditions. 

Hence (X J) . _ denotes X at reference conditions with the exception of 

ref A, B 

A and B which are non-reference. It follows that the change in payload ratio 
from reference conditions due to a change in A and B is 

AX A, B = 'XefU, B ■ "ref = * I AA(F/W o> + AB f^)] (10) 


Changes in X due to variations in T , for a given I , can 
, , tr sp 

be expressed as 


(AX)_ = (X 4 - X r 

v 'T v ref'T ref 

tr tr 


( 11 ) 
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Calculations were made varying T^. in increments of 5, 10, and 20 hours 

over the range considered. The changes in payload ratio were computed 
using Equation 11 and are shown in Figure 9. Linear interpolation between 
Tt r curves is sufficiently accurate for preliminary design purposes. 

The values of are independent of the vehicle specific weight 

tr 

parameters assumed, but are slightly dependent upon initial orbit altitude. 
However, it was found that within the altitude range considered, this 
dependency can be neglected. 


Changes in the initial orbit altitude, other reference conditions 
remaining constant, also must be accounted for in estimating payload ratio. 

It was discovered that the change in payload ratio for an incremental change in 
altitude is approximately the same over the range of altitude considered, 
regardless of the value of I . Consequently, altitude and transfer time 

corrections to the reference payload ratio can be made independently. 

The same applies to the other parameters. 


Generally initial orbit altitudes will be lower than 555. 6 km 

(300 N. M. ); therefore the corrections for initial altitudes will be based 

upon a 185. 2-km (100-N. M. ) incremental decrease. The change in 

reference payload ratio for an incremental decrease of 185. 2 km (100 N. M. ) 

is shown in Figure 10. Linear interpolation between curves presented 

for the assumed value of I is sufficient, since altitude corrections are 

sp 

generally small. For arbitrary initial altitude, h , the change in payload 
is given by 


(AM h 

O 


(h ) , 

o ref 


185. 



555. 6 -h 

(AX h * 185.2 = 185. 2 

o 


where (AX, ) is the correction shown in Figure 
V h ' 185.2 6 

o 

are measured in km. 


(AX ) 

h 185.2 
O 


10, (h ) 

o ref 


( 12 ) 


and h 

o 
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Combining the corrections given by Equations 10, 11, and 
12, the payload ratio for any assumed set of conditions can be found by 


A, B, h_ = X__ r + AX a ^ + AX,p + AX^ 


tr o ref A, B 


tr 


The effect of non-zero values of C can now be shown. Since 
C is a function of W only, it can be subtracted directly from (X) 


C = 0‘ 


Hence 


X A, B, C, T , h X A, B, T , h ' C 
tr o tr o 


° r X A, B,C,T , h X ref + AX A, B + AX T + A \ 
tr o tr o 


(13) 


Using this relation, the payload ratio for a particular vehicle with a 
given mission can be computed. 


B. CHARACTERISTIC VELOCITY 

The characteristic velocity, AV.^, representing a 

combination of the inertial velocity requirement and gravity losses, is 
computed in a manner simular to that described for estimating payload 
ratio. Reference curves for AV\^ are presented in Figure 3. Correction 

curves for change in transfer time requirement and initial orbit altitude 
are presented in Figure 11 and 12, respectively. The corrections for 
Tj. r were computed from the relation, 


A < AV id>T = I (AV id>refl T ‘ < AV id>ref 
tr tr 

which is analogous to Equation 11. 


(14) 
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The change in reference characteristic velocity for an 
incremental altitude decrease of 185. 2 km (100 N. M. ) is shown in 
Figure 11 and corrections for arbitrary initial altitudes are given by. 


(h ) -h 

(AV.,), = -- iaV 7"° [A(AV.,), ] 

id h 185.2 L id h 185.2 

o o 


555. 6 - h 

■ ■■ - ° [A(AV. ,), ] 

185. 2 id h 185.2 


(15) 


with [ A(AV.^ ] 2 determined from Figure 12. 


given by 


The characteristic velocity for arbitrary T and h is then 

tr o 


< AV id>T ,h = < AV id>ref + A < AV id>x +A(AV id)h 
tr o tr o 


( 16 ) 


C. PROPELLANT RATIO 

With the value of AV. , corrected to the desired conditions, 

id 

the propellant ratio may be determined using the familiar relation 


A V. = g I In - 

id ®n sp 1 - C 


rearranging, 


(U T ,h = 1 ' eXP 
tr o 


- (AV id>T , h 
tr o 

g I 
6 n sp 


(17) 


The reference curves for shown in Figure 4, are used 
to determine payload ratio, in Equation 10. 
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D. BURNING TIME 


As in the calculation of burning time is readily determined 
from Equation 6. Hence, 


I <&) T 

sp T 


, h 

tr o 


. , h 
tr o 


(F/W ) 
o 


(18) 


The reference burning time curves in Figure 5 are shown 
only for the purpose of illustrating the typical variation with F/W Q 
and I 

sp 


E. FLIGHT PATH ANGLE AND CENTRAL ANGLE 

The procedures for estimating the flight path angle and the 

central angle are identical to the procedure outlined for calculating 

AV. , and they are discussed together. Thus, analogous to Equations 14 

and 15, the corrections for 0 and ip with T are given by 
' t r 


(A 6 ) 

tr 


^ref^T 


- 0 


tr 


ref 


and 


(a^) t 

tr 


^ref^T 


- 4^ 


tr 


ref 


(19) 


( 20 ) 


These corrections are shown in Figures 13 and 15, 
re spectively. 

The corrections for 6 and ^ for arbitrary h are determined 
from the relations, 


(a e\ 


555. 6 - h 

185. 2 * A0 h ^ 185.2 
o 


and 


(A^h 


555. 6 - h 

185. 2 * 185.2 

o 


( 21 ) 


( 22 ) 
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where (A 9 , ) and (Aik ) > the corrections for a 185.2-km (100-N. M. ) 

h 185.2 h 185.2 

o o 

incremental decrease in initial orbit altitude, are given in Figures 14 and 16. 
It was discovered that the initial orbit altitude corrections for these 
parameters were insensitive to specific impulse variations. Hence, the 
resulting corrections are presented by only one curve for each parameter. 

It follows that, for arbitrary conditions, 


,h =0 ref +(A(?) T +(A6) h 

tr o tr o 


and 


(40 , = 4 r + (A4)m + A4k 

V T , h T ref v T v 'h 
tr o tr o 


(23) 


(24) 


F. ALTITUDE AT INJECTION 


The procedure for estimating the altitude at injection is the 
same as for the previous parameters. Thus, the correction for T is 
given by 


(A 1 Vt ^b^ref^T ' ^b^ref 
tr tr 


where (Ah ) is shown in Figure 17. 
b T 

tr 


(25) 


The correction for h for arbitrary h is determined from 
. . b o 

the relation, 


< A Vh 

< 

where [ (Ah^)^ 


555. 6 - h 

= 185! 2 ” ^ Ah bV ^ 185.2 

o 

1 „ , the correction for 185.2-km incremental decrease 

J 185.2 ’ 


o 

of initial altitude, is given in Figure 18. 
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Then for arbitrary conditions, 

( Vt ,h = (h, ) + (Ah. )_ + (Ah ) (27) 

tr o b ref b T b h. 

tr o 

This completes the estimation procedures for all the parameters 

considered. 


SECTION IV. ACCURACY EVALUATION 

Results of several numerical examples are compared with exact 
computer solutions and the associated errors for each parameter are 
included in Table 2. Two of the numerical examples are presented in 
detail in Appendix A for the purpose of illustrating the performance 
estimation procedure. One example is included, in which the initial 
orbit altitude is outside the range assumed, to determine the accuracy 
with which the data can be extrapolated. 

The assumed data for the numerical examples were deliberately 
chosen so that interpolation was necessary in most cases. Hence, the 
results should be indicative of the maximum errors which can be expected 
in practical applications of the procedure. Since the examples are not 
numerous enough to represent a statistical sampling, the errors shown 
should not be construed as being the most probable. 


SECTION V. CONCLUSIONS 

1. Several numerical examples were worked to test the procedure 

for both speed and accuracy. From the experience gained by working these 
examples it is concluded that the procedure is rapid and has sufficient 
accuracy for preliminary design studies. 

2. From the results of the numerical examples and their associated 
accuracy, presented in Table 2, it is concluded that the average errors 
that can be expected in a particular application will be less -than 1 per cent 
for all the parameters considered. The maximum errors associated with 
the various parameters that can be expected with careful use of the procedure 
are as follows: 
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< 1 % 


AV id> c, t b , e, qj 




< 2 % 


The burnout altitude estimation exceeds 1 per cent only for the cases 
involving low thrust -to -weight ratios; however, the accuracy is still 
within preliminary design requirements. 


3. Performance estimations can be made for any nuclear vehicle 
and any lunar mission profile, within the range assumed for the associated 
parameters. Also, this method of performance evaluation should have the 
capability of being extended to any orbit-launched systems with equivalent 
success. 
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APPENDIX A 


NUMERICAL EXAMPLES 


The purpose of these numerical examples is to illustrate the 
methods utilized in the estimation procedure. Five numerical examples 
were solved, and the results are presented in Table 2; however, only- 
two of these are presented here. 

A. NUMERICAL EXAMPLE 1 


The first step in'the estimation procedure is the determination 

of reference values, since these are the basis for the procedure. The 

arbitrarily assumed conditions for Example 1 are presented below. At 

the values of F/W^ -0.2 and I = 735 seconds, the necessary reference 
° sp 

values can be read from Figures 2, 3, 4, 6, and 7 as follows: 

As sumed Conditions Reference Values 


F/W q = 0.2 

I = 735 s 
sp 

A = 0. 2 

B = 0. 15 

C = 0. 0 

T = 85 hr 
tr 

h = 277. 8 km 


X , = 0. 5863 
ref 

C , = 0. 3579 
ref 

(AV. ,) r = 3197 m/s 
id ref 

^ = 1. 5275 rad (87. 5°) 

ref 

6 = 1. 1572 rad '{66. 3°) 

ref 

(h ) f = 2 300 km 
b ref 


The first performance parameter discussed, \, is probably 
the most important one. Corrections of \ ^ are necessary for changes 

in A, B, C, T , and h . For the assumed values of A and B, 
tr o 

Equations 8 and 9 yield: 


A A = A - 0. 10 


= 0 . 20 - 0 . 10 


= + 0 . 10 


15 


AB = E - 0. 10 


= 0. 15 - 0. 10 
= +0.05 

From Equation 10 the change in \ ^ due to non-reference 

value s of A and Bis 

AX A,B = " l* A (F/W o ) + AB (£ ref )] 

= - [0.10 (0. 20) + 0. 05 (0. 3579)] 

= -0. 0379 

From Figure 9 the change in K _ due to non-reference T 

ref tr 

is 

A\ t = +0. 0017 


From Figure 10 and Equation 12 the change due to non- 
reference h is 
o 


(A\_ ) = -0. 0056 

h 185.2 

o 


from Equation 12 


A \ 


o 


555. 6 - h 

185.2 (AX h \ 8 5.2 
o 


A \ 


o 


555. 6 - 277. 8 
185. 2 


(-0. 0056) 


= -0. 0084 
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Therefore the corrected value of \ is 

X = X ref + AX A, B + A V +AX h " C 

tr o 

= 0.5863 - 0.0379 + 0.0017 - 0.0084 - 0 
= +0. 5417 

The characteristic velocity is determined in a similar 

manner. The T correction is found in Figure 11 to be 
tr 


A( Av id ) T = -16 m/s 


tr 


From Figure 12 and Equation 15 the correction in (AV. d ) rg{ 


for non-reference h is 

o 

[ A(AV. ,), ] = + 55. 8 m/s 

L id h j 185.2 

o 


from Equation 15 
A(AV. ,) 


555. 6 - h 

id'h = 185.2 t A(AV idV ^ 185.2 

o O 


555. 6 - 277. 8 


(55. 8) 


185. 2 

= +83. 7 m/ s 

Therefore from Equation 16, 

AV. j = (AV. d ) ref + A(AV ; „) T + A(AV. d ) 


id 


id T 


tr 


id h 


= 3197 - 16 + 83. 7 
= +3265 m/ s 
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The propellant ratio is determined from Equation 17 as 


follows, 


t, = 1 - exp 


-< AV id>T ,h 
tr o 

x sp 


= 1 - exp 
= +0. 3638 


-3265 


9. 82 (735) 




The burning time is found from Equation 18, 

1 (4) T h 

SP T tr* h o 
F/W 


735 (0. 3638) 

0 . 20 


= +1337 sec 


The T correction for 6 .is found in Figure 13 to be, 
tr ref 


A Orj, = +0. 0019 rad 

tr 

From Figure 14 and Equation 21 the change in 0 ^ i- s » 


(A0 h \s5.z = -°- 0158 rad 


and from Equation 21, 

555.6 - h 

A 9, = fnr '- T ~ ~ (A 6. ) 

h 185.2 h 185.2 

o o 


555. 6 - 277. 8 
185.2 

= -0. 0236 rad 


(-0. 0158) 
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The corrected value of 6 , from Equation 23, is then 

6=6 . + A 6 + A 6 

ref T h 

tr o 

= 1. 5272 + 0. 0019 - 0. 0236 
= +1. 1355 rad (65.06°) 

The central angle at injection is determined by the same 

method. The T correction read from Figure 15 is: 
tr 


A'I'rp = -0.0063 rad 
tr 


From Figure 16 and Equation 22 the change in 4* ^ for 

non-reference h is: 
o 

(A^V ) = +0. 0705 rad 

h 185.2 
o 


and from Equation 22, 

555. 6 - h 

A4V = 


o 


185. 2 ^185.2 

o 


555.6 - 277. 8 
185. 2 


(+0. 0705) 


= +0 . 1058 rad 

The corrected value of from Equation 24, is 

+ =+ ref +A+ T +A* h 

tr o 

= 1. 5272 - 0. 0063 + 0. 1058 
= +1. 6267 rad (93. 20°) 
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The burnout altitude is determined in a similar way. The 
burnout altitude correction for non-reference transfer time is found in 
Figure 17 to be 

(Ah b ) T = -21 km 
tr 

The correction for the non-reference h is determined by 
using Figure 18 and Equation 26. From Figure 18 

[ (Ah ) ] = -68 km 

b h 185-2 

o 


and from Equation 26, 

555 . 6 - h 


^ Ah bV 185.2 

o 


^ Ah b*h ] 185.2 


555. 6 - 277. 8 
185. 2 


(- 68 ) 


= -102 km 


Therefore from Equation 2 7 the corrected value of h^ is 


\ (1 Vref + (A Vt + (Al Vh 

tr ( 

= 2300 - 21 - 102 


= +2177 km 


B. NUMERICAL EXAMPLE 2 

The arbitrarily assumed conditions for Example 2 are presented 
below. At the values of F/W Q = 0.4 and I S p = 980 seconds, the necessary 
reference values can be read from Figures 2, 3, 4, 6, and 7 as follows: 
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Assumed Conditions 


Reference Values 


F/W q = 0.4 

I = 980 s 
sp 

A = 0. 09 

B = 0. 15 

C = 0. 02 

T = 78 hr 
tr 

h = 463.0 km 


\ = 0. 6565 

ref 

^ e £=°- 2756 

(AV id ) ref= 3107 m/s 
4* f = 0- 8465 rad (48. 5°) 

9 r = 1. 3308 rad (76. 2°) 
ref 

(h ) = 1090 km 

b ref 


Corrections for X are necessary for changes in A, B, C, 

T tr' and h o' F ° r the assume d values of A and B, Equations 8 and 9 
yield: 


AA = A - 0. 10 
= 0. 09 - 0 . 10 
= - 0.01 
AB = B - 0. 10 
= 0.15- 0. 10 
= +0. 05 

From Equation 10 the change due to A and B is 
AX A,B = -t AA < F/W o> + AS (C ref )] 

= - [ -0. 01 (0. 4) + 0. 05 (0. 2756)] 

= - 0. 0098 
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From Figure 9 the T correction is 

tr 

A\ t = +0. 0008 
tr 


From Figure 10 and Equation 12 the h^ correction is 
(AX, ) = -0. 0041 

h 185.2 
O 


and from Equation 12, 

555. 6 - h 
AX, = — , -TT A -^ 


185.2 * AX h \ 8 5.; 


555. 6 - 463. 0 
185. 2 


(-0. 0041) 


= - 0 . 0021 


Therefore, the corrected value of X. is 


X = X , + AX a + AX + AX -C 
ref A, B T h 

tr o 


= 0. 6565 - 0. 0098 + 0. 0008 - 0. 0021 - 0. 0200 


= +0. 6254 


The characteristic velocity is determined in a similar 
manner. From Figure 11 the T correction is 

A(AV id ) T = _ 8> g m / s 

tr 

From Figure 12 and Equation 15 the h Q correction is 


[A(AV.j), ] = +45.5 m/ s 

L id h J 185.2 
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and from Equation 15, 

555. 6 - h 


A(AV iA 


185.2 ^ A(AV id^h ^ 185 .2 


555.6 - 463 . 0 
185. "2 


(45. 5) 


= +22.8 m / s 


Therefore from Equation 16 the corrected value of AV.^ 

AV id = (AV id ) ref + A(AV id>T + A(AV id ) h 

tr o 

= 3107 - 8. 8 + 22. 8 
= +3 121 m / s 


The propellant ratio is determined from Equation 17 as 


follow s : 


t, = 1 - exp 


= 1 - exp 


(AV id^T , h 
tr o 


g I 

B n sp 


-312 1 

(9. 82) (980) 


= +0. 2775 


t b = 


From Equation 18 the burning time is 

VA ,h 

tr o 


F/W 


980 (0. 2775) 
0.4 


= +679. 9 sec 



IS 


In Figure 13 the T correction for 0 is 

tr 

A6 t =+0.0061 rad 
tr 

From Figure 14 and Equation 21 the change in 0 
(A0, ) = -0.0105 rad 

h 185.2 

o 

and from Equation 21, 


A & 


555. 6 - h 

( 

185. 2 “ 


(AO, ) 

h 185.2 
O 


555. 6 - 463. 0 
185. 2 


(-0. 0105) 


= -0. 0052 rad 

The corrected value of 6 is then, 


9 = 6 . + A0_ + A0, 

ref T h 

tr o 


= 1. 3328 + 0. 0061 - 0. 0052 


= +1. 3263 rad (76. 0°) 

The central angle at injection is determined by a similar 

method. From Figure 15 the T correction is 

tr 

A^rp = -0. 0028 rad 
T tr 

From Figure 16 and Equation 22 the h^ correction is 

< A \ W* = +0-0422 rad 
o 
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and from Equation 22, 

555. 6 - h 

Aik = 


185.2 *185.2 


o o 

555. 6 - 463. 0 


185. 2 


(0. 0422) 


= +0. 0211 rad 

The corrected value of 4 1 is then, 

^ ^ , + A 4* + A 4*, 

ref T h 

tr o 

= 0. 8465 - 0. 0028 + 0. 0211 
= +0. 8648 rad (49. 55°) 

The burnout altitude is determined in a similar manner 
The transfer time correction from Figure 17 is 

(Ah b ) T =-5.0 km 
tr 


The correction for arbitrary initial altitude, from Figure 18 
and Equation 26 is 

(Ah ) = -144 km 

b h 185.2 

o 

and from Equation 26 

555. 6 - h 

< A Vh = 185.2 ° t (A Vh 1 

o o 


555. 6 - 463. 0 
185. 2 


(-144) 


= - 7 2 km 
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The corrected value of from Equation 2 7 is then 


h b ^b'ref + 


< a Vt + (A Vh 

tr o 


= 1090 - 5. 0 - 72 
= +1013 km 


The results of these examples, along with three others are 
compared with exact computer runs and the per cent errors associated 
with each of the parameters are presented in Table 2. 


26 



Inertial Arrival Velocity At Periselenum, km/ s 



- HRS FROM INJECTION 




Reference Payload Ratio, 


0. 6800 
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FIGURE 2. REFERENCE PAYLOAD RATIO 




Thrust -to -Weight Ratio, F/W Q 


FIGURE 3. REFERENCE CHARACTERISTIC VELOCITY 






Reference Propellant Ratio, 
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FIGURE 4. REFERENCE PROPELLANT RATIO 
30 
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FIGURE 5. REFERENCE BURNING TIME TO INJECTION 
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FIGURE 7. REFERENCE CENTRAL ANGLE AT INJECTION 
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Reference Burnout Altitude, (h ) 
















Change In Reference Payload Ratio, AX., 
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FIGURE 9. CHANGE IN REFERENCE PAYLOAD RATIO FOR NON-REFERENCE 
TRANSFER TIMES 
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FIGURE 10. CHANGE IN REFERENCE PAYLOAD RATIO FOR INCREMENTAL ALTITUDE 

DECREASE OF 185. 2 KM ( 100 N. M. ) 






Change In Reference Characteristic 
Velocity, [A(AV. d ) h 1 , m/s 
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FIGURE 12. CHANGE IN REFERENCE CHARACTERISTIC VELOCITY 
FOR INCREMENTAL ALTITUDE DECREASE OF 185. 2 KM (100 N. M. ) 
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Change In Reference Flight Path Angle, AO, 


h"* ( rad ) 
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13. CHANGE IN REFERENCE FLIGHT PATH ANGLE FOR NON-REFERENCE TRANSFER TIMES 
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Change In Reference Central Angle, A4f 
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FIGURE 15. CHANGE IN REFERENCE CENTRAL ANGLE FOR NON-REFERENCE TRANSFER TIMES 
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FIGURE 16. 


CHANGE IN REFERENCE CENTRAL ANGLE FOR INCREMENTAL ALTITUDE DECREASE 
OF 185. 2 KM (100 N. M. ) 
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Change In Reference Burnout Altitude, (Ah ) 




Table 1. Data Assumed For Numerical Examples 


Example 

Number 

1 

2 

3 

4 

5 

1 (sec) 

sp 

735 

980 

750 

755 

885 

F/W 

o 

0. 20 

0. 40 

0. 30 

0.45 

0. 10 

h (km) 
o 

277. 8 

463. 0 

185. 2 

835. 0 

694. 9 

T tr< hr > 

85 

78 

110 

57 

60 

A 

0. 20 

0. 09 

0. 05 

0.12 

0. 08 

B 

0. 15 

0. 15 

0. 12 

0. 22 

0. 07 

C 

0. 00 

0. 02 

0. 03 

0. 00 

0. 00 
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Table 2. Comparison of Exact and Predicted Performance Parameters 
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